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Abstract 

We report the results of a very high statistics Monte Carlo study of 
the continuum limit of the two dimensional 0(3) non-linear a model. 
We find a significant discrepancy between the continuum extrapolation 
of our data and the form factor prediction of Balog and Niedermaier, 
inspired by the Zamolodchikovs' S'-matrix ansatz. On the other hand 
our results for the 0(3) and the dodecahedron model are consistent 
with our earlier finding that the two models possess the same contin- 
uum limit. 



In a recent paper [||] we reported some striking numerical results: the 
continuum limit of the two dimensional (2D) 0(3) nonlinear a model seems 
to agree as well with the form factor prediction inpired by Zamolodchikovs' 
S'-matrix as with the continuum limit of the dodecahedron spin model. The 
latter, known rigorously to possess a phase transition at nonzero temperature, 
is almost certainly not asymptotically free. Zamolodchikovs' ansatz on the 
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other hand is beheved to embody asymptotic freedom, a behef which was 
reenforced by some remarkable scahng properties discovered by Balog and 
Niedermaier Q]. 

Initially ||^ we thought that the resolution of this apparent paradox must 
come from the absence of asymptotic freedom in the form factor approach 
and that the scaling proposed by Balog and Niedermaier must actually be 
false. Soon after though we realized that the trouble was much more serious: 
indeed from the direct computation of Balog and Niedermaier of the 2,4 and 
6 particle contribution it followed that the asymptotic value of the transverse 
current two point function J{p) must obey the inequality 



0. On the other hand in two recent papers |^, || we used reflection positivity 
to prove rigorously an upper bound on J(oo) in terms of Pert- In those papers 
the bound is given only for the 0(2) model, but it is straightforward to 
generalize it to any 0{N). For 0(3) it reads 



Combining these two inequalities one would conclude that Pert > 5.5. This 
lower bound is so large that irrespective of the Balog and Niedermaier scaling 
ansatz, it suggests that the critical point may very well be at infinity and the 
model asymptotically free. 

Since, as stated above, asymptotic freedom of 0(3) is hard to reconcile 
with its having the same continuum limit as the dodecahedron spin model, 
we decided that we needed a more refined numerical study to decide these 
issues. In the present paper we will report the results of our investigation 
regarding the continuum limit of the 0(3) and dodecahedron spin models. 
We strove to achieve very high numerical accuracy and to our knowledge our 
accuracy, exceeding sometimes .05%, has never been achieved before. Our 
results suggest that in fact the model defined by the form factor approach 
does not agree with the continuum limit of the 0(3) model. On the other 
hand, the numerical evidence is consistent with the hypothesis that the 0(3) 
and the dodecahedron spin models have the same continuum limit. 

We begin our report by specifying the models and defining the quantities 
measured. In both models the action is nearest neighbour, the lattice square 
of length L, with periodic boundary conditions and the inverse temperature 
p. The quantities measured are (with P = (p, 0)): 



J(oo) > 3.651 



(1) 



J{00) < -Pert 



(2) 
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• Spin 2-point function 

G{p) ^ ^,mm; ^(pj^Ee^^'M^) (3) 

• Current 2-point function J(p) = F(0) — F{p) with 

Fip) ^ ^,j:mpn, j;ip)-i:e'^xi^) (4) 

a=l X 

where 

f{x) = l3eabcSb{x)sc{x + /i) 



• Magnetic susceptibihty 

X = G'(O) (5) 

• Effective correlation length ^: 

1 



2sin(7r/L) 



(G(0)/G(1) - 1) (6) 



Zero momentum renormalized coupling: 



where is the invariant connected spin 4-point function at zero mo- 
mentum. 

Some comments regarding these quantities: 

• Our definition of ^ differs from that of Balog and Niedermaier, who 
define 1 /^^ as the position of the pole in the spin 2-point function (on 

our Euclidean lattice, that pole is inaccessible). If wc implement our 
definition upon the form factor prediction, using = 1 as the lowest 
nonvanishing momentum, we find ^ = 1.0016187, compared to 1. 

• The renormalized spin 2-point function is Gr(pO — ^(pO/^(0) (nor- 
malized so that Gr(0) — 1). 
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• For the 0(3) model the definition of the current, including the factor 
of /3, is fixed by the Ward identity. The quantity Jo{3){pO is a 'renor- 
malization group invariant', which describes continuum physics in the 
limit ^ — oo, fixed. 

• For the dodecahedron spin model the definition of the current is ad 
hoc, but if in fact its continuum limit agrees with that of 0(3) it is 
possible that upon multiphcation by a suitable factor c(^), c{^)Jdod{pO 
reaches the same limit as Jo{3){pO- The results reported in [jl| seemed 
to corroborate this assumption. 

• Besides being a renormalization group invariant, Qr vanishes in any free 
field theory. Hence Qr tests the nontrivial part of the S'-matrix. 

The numerics were performed on a SGI-2000 machine. Both models were 
simulated using the cluster algorithm and all quantities, except gr, measured 
using Wolff's improved estimator |Q. At every jS and L value we performed 
runs consisting of 100,000 thermalization clusters and 1,000,000 clusters used 
for measurements. These runs were repeated for 0(3) up to 149 times and for 
the dodecahedron up to 246 times. We give the results for x and ^ in Tab.l 
(0(3)) and Tab. 2 (dodecahedron). As can be seen, our results for 0(3) are 
more precise than for the dodecahedron, in spite of the extra effort invested 
in the latter. The errors were estimated using the jack-knife method on the 
results of different runs, since this led to a larger estimated error than the 
error estimated in each run. 

The form factor prediction gives the continuum, infinite volume values 
of Gr{p^) and J{pC,) (however truncated as to the number of intermediate 
particles included). To relate Monte Carlo measurements to these predictions 
one must take both the continuum limit, ^ — >■ cx3, and the thermodynamic 
limit, L/^ — > oo. In Fig.|I] we present the deviations from the continuum, 
infinite volume limit for the free field, displaying J(p^) for several correlation 
lengths ^ and ratios L/C,- It can be seen that for pC, < 12 for L/C, > 13 the 
finite volume corrections become negligible (less than 2 x 10~^) for ^ > 11 
(the smallest correlation length in our study). We verified directly that the 
same is true for the 0(3) model by running at the same ^ values on L/C 
approximately 7-8, 13-15 and 23-25. 

Another lesson learned from Fig.|I] concerns the approach to the contin- 
uum: at fixed pC the approach is described by 

a/e+b\ogic)/e (8) 
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with a and b depending upon p^. Since at small p$, the 0(3) data are rea- 
sonably close to the free field values, we employed this ansatz in all our 
extrapolations to the continuum limit. 




Figure 1: Lattice and finite volume artefacts for the complex free field; dotted 
curves: L/^ ^7, solid lines: L/^ ^ 14, dashed lines: L/^ ^ 24 - the latter 
two are practically indistinguishable. 

As stated above, our definition of ^ differs from that of Balog and Nie- 
dermaier. Even though we could correct for this difference by applying our 
definition to their data, we decided that a better approach is to eliminate ^ 
altogether and plot J versus G{0)/G{p). For small p^ the latter quantity is 
very close to its free field value, 1 + {p^Y- This plot is shown in Fig.H for 
various ^ values and L/^ aproximately 14. The solid lines represent a least 
square fit to the data of the form 

Ea,(J;,ee(K))' (9) 
k=l 

where 



5 



0.06 




0.04 



0.02 



__L^140,(=1 1.06 



--i 



J,=250,(=19.00 
- t — 



--t --f 



5 10 
G(0)/G(p)-1 



15 



Figure 2: Comparison of Monte-Carlo (Jo(3)) with form factor results (J//); 
the dashed line represents the continuum extrapolation. 



with z = (pO^ is the expression corresponding to J in the free two-component 
scalar field theory. We also show the extrapolation of our results to the 
continuum limit using the ansatz in eq.(3) and the form factor prediction. 
As it can be seen, the latter disagrees with the former. 

The disagreement is statistically significant, yet very small. One may 
wonder how such a wonderful agreement could occur if the form factor pre- 
diction is actually false. A possible answer is this: the correlation functions 
J and Gr are, especially at low momenta, very close to those of the free 
theory, deviating only by less than 5%. The reason behind this is that the 
form factor squares contributing to the spectral functions of the correlation 
functions are at low momenta 

• dominated by the lowest possible number of intermediate particles (1 
for Gr, 2 for J) 

• the squares of those lowest form factors are very close (in the case of 
Gr even equal) to those of the free theory. 

Therefore, if we want to test the form factor prediction, this universal free 
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Figure 3: Comparison of Monte-Carlo ( Jo(3) ^ dotted lines) and their contin- 
uum extrapolation (dashed line) with form factor results {Jjf - solid line): 
interacting contribution. 

field value should be subtracted. To achieve such a subtraction we do the 
following: we plot x'^J{x) versus G{0)/G{x) x = for both the 0(3) model 
and the free field. We then subtract the two curves at the same value of the 
abscissa and plot the difference. This is shown in Fig.^, together with the 
form factor prediction. The latter differs from the continuum value extrapo- 
lation by about 25% at p = 0. 

Another quantity which measures directly the nontrivial part of the 5*- 
matrix is the renormalized coupling Qr- For the 0(3) model this quantity was 
determined to be around 6.5-6.9 and we determined it again at /3 = 1.5, 
^ = 11.06 to gr = 6.565(25) and at /? = 1.6 to gr = 6.609(18). Since gr is 
expected to be monotonically decreasing in j3, the small difference between 
the two values is probably just due to statistical fluctuations. The value of 
gr in the form factor approach has not yet been calculated, even though the 
ingredients are the same as for G{p). We think that the comparison would 
be most interesting. 

Next we would like to present the comparison of the dodecahedron spin 
model with the 0(3) model. If the two models possess the same continuum 




Figure 4: Comparison of Gj. for 0(3) and dodecahedron. Data with ^ ^ 
11, 19, 34, 65 are overlaid. The sohd hnes are fits. 

hmit, then one would expect that for sufficiently large ^ both the lattice 
artefacts and the finite volume corrections become identical. Consequently 
if we use the same definition for ^ in both models, it should be legitimate to 
compare the same 'renormalization group invariants' at the same pC, value. 
In Fig.| we present the comparison of Gr- The data are consistent with the 
hypothesis that the two models have the same continuum limit. The error 
bars are dominated by the those coming from the dodecahedron, where the 
fiuctuations are very large. We present a similar comparison in Fig.^ for J, 
where the renormalization was chosen so that the rms deviation between the 
two models in the range from = to = 15 is minimized. As stated 
above, the theoretical status of this procedure is not clear. Finally we report 
on gr for the dodecahdron. We measured it at ^ = 11 and 19 and obtained 
6.42(84) and 6.33(83). These values suggest that its continuum limit value 
is also around 6.5, in agreement with the value of Qr for 0(3). 

Our conclusion is that whereas there is good numerical evidence that 
the model defined via the form factor approach does not agree with the 
continuum limit of the lattice 0(3) model, there is no evidence that this 
limit is different from that of the dodecahedron model. These facts represent 
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Figure 5: Comparison of J for 0(3) and dodecahedron (renormalized) . Data 
with ^ ^ 11, 19, 34, 65 are overlaid. The sohd hnes are fits. 

a major setback for the accepted saga regarding the special properties of the 
nonlinear a models with > 3. 

We acknowledge numerous discussions with J.Balog and M.Niedermaier. 
A. P. is grateful to the Max-Planck- Institut for its hospitality. 

Tab.l: 



Susceptibility x ond effective correlation length C, for 0(3) 



p 


L 


# of runs 




X 


1.5 


140 


149 


11.0600(29) 


176.7348(566) 


1.6 


250 


86 


18.9971(55) 


448.3651(1515) 


1.7 


500 


100 


34.5712(97) 


1270.8185(4196) 


1.8 


910 


44 


64.8723(277) 


3839.8609(2.1498) 



Tab.2: 

Susceptibility x o,nd effective correlation length ^ for the dodecahedron 
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L 


# of runs 




X 


1.49 


140 


246 


11.0540(40) 


176.7967(840) 


1.5835 


250 


201 


19.0337(55) 


451.1633(2320) 


1.672 


500 


139 


34.4139(167) 


1265.2500(7350) 


1.76 


910 


78 


66.6195(522) 


4046.1542(3.9885) 



References 

[1] A.Patrascioiu and E.Seiler, Phys. Lett. B 430 (1998) 314. 

[2] J.Balog and M.Niedermaier, Phys. Lett. B 386 (1996) 224. 

[3] J.Balog and M.Niedermaier, Nucl Phys. B 500 (1997) 421. 

[4] J.Balog and M.Niedermaier, Phys. Rev. Lett. 78 (1997) 4151. 

[5] A.Patrascioiu and E.Seiler, Phys. Rev. Lett. 80 (1998) 639. 

[6] U.Wolff, Phys.Rev.Lett. 62 (1989) 361. 

[7] A.Patrascioiu and E.Seiler, Phys.Lett.B 417 (1998) 123. 

[8] A.Patrascioiu and E.Seiler, Phys.Rev. E 57 (1998) 111. 

[9] A.Patrascioiu and E.Seiler, Nuovo Cim. 107A (1994) 765. 



10 



